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A consistent comparison between nonrelativistic and relativistic descriptions of the (e, e'p) re- 
action is presented. We use the nonrelativistic DWEEPY code and develop a fully relativistic model 
starting from the effective Pauli reduction for the scattering state and the relativistic mean field the- 
ory for the bound state. Results for the 0(e, e'p) differential cross section and structure functions 
are compared in various kinematical conditions. A limit in energy of the validity of the nonrelativistic 
approach is established. The effects of spinor distortion and of the effective momentum approxima- 
tion for the scattering state are discussed. A satisfactory agreement with data of differential cross 
sections, structure functions, and polarization observables is achieved. 

PACS numbers: 25.30.Fj, 24.fO.Jv, 24.70. +s. 



I. INTRODUCTION 



Exclusive (e, e'p) knockout reactions represent a very clean tool to explore the single-particle (s.p.) aspects of 
nuclear structure revealing the properties of proton- hole states contained in the hole spectral function (TJ-Q] . 
^ ' Several high-resolution experiments were carried out at Saclay J|§ and NIKHEF §]. The analysis of the missing 
energy and momentum dependence of the experimental cross sections allowed to assign specific quantum numbers and 
spectroscopic factors to the peaks in the energy spectrum. The calculations for this analysis were carried out with 
the program DWEEPY 0, within the theoretical framework of a nonrelativistic distorted wave impulse approximation 
(DWIA), where final-state interactions (FSI) and Coulomb distortion of the electron wave functions are taken into 
y—\ ', account. Phenomenological ingredients were used to compute bound and scattering states. The outgoing nucleon 
• scattering wave functions are eigenfunctions of an optical potential determined through a fit to elastic nucleon- 
nucleus scattering data including cross sections and polarizations. Bound-state wave functions were calculated with a 
Woods-Saxon well, where the radius was determined to fit the experimental momentum distributions and the depth 
was adjusted to give the experimentally observed separation energy of the bound final state. This theoretical approach 
was able to describe, with a high degree of accuracy, in a wide range of nuclei and in different kinematics, the shape 
of the experimental momentum distributions at missing-energy values corresponding to specific peaks in the energy 
spectrum. In order to reproduce the size of the experimental cross sections, the normalization of the bound-state 
wave function was fitted to the data and identified with the spectroscopic factor. These values, however, are smaller 
than those predicted by many-body theories. 

Similar models based on a fully relativistic DWIA (RDWIA) framework were developed in more recent years p~4| . 
In these approaches the bound nucleons are described by s.p. Dirac wave functions in the presence of scalar and 
vector potentials fitted to the ground-state properties of the nucleus, and the scattering wave function is solution 
of the Dirac equation with relativistic optical potentials obtained by fitting elastic proton-nucleus scattering data. 
Some of these approaches include also an exact treatment of the Coulomb distortion of the electron waves P,[l0|]. 
Also RDWIA analyses were able to give a good description of the experimental momentum distributions. Thus, the 
shape of the distributions in DWIA and RDWIA calculations are similar, while the spectroscopic factors, obtained by 
scaling the calculated cross sections to the data, are in RDWIA about 10 — 20% larger than in DWIA analyses, and 
thus closer to theoretical predictions. The difference was attributed to the relativistic optical potential, which leads to 
stronger absorption. The nucleon-nucleus interaction exhibits characteristic non-localities which arise quite naturally 
in the Dirac approach and whose effect was not included in standard nonrelativistic DWIA calculations, but can be 
accounted for in the nonrelativistic treatment by a renormalization of the scattering wave function. The so-called 
Darwin normalization factor |f~5| l , that essentially changes the Schrodinger wave function into the upper component of 
the Dirac wave function, increases the absorption due to FSI and produces a quenching of the calculated cross section 
by about 15%, with a corresponding enhancement of the spectroscopic factor in agreement with the results obtained 
in RDWIA. 

New data have recently become available from TJNAF. The cross section for quasielastic fp-shell proton knockout 
has been measured and the response functions and the asymmetry have been extracted in the 16 0(e,e'p) reaction at 
four-momentum transfer squared Q 2 — 0.8 (GeV/c) 2 and energy transfer u ~ 439 MeV ]lr| . In the same kinematics 



1 



also first polarization transfer measurements have been carried out for the 16 0(e, e'p) reaction u7\- The polarization of 
the ejected proton in the 12 C(e, e'p) reaction has been measured at Bates with Q 2 = 0.5 (GeV/cp and outgoing-proton 
energy T p = 27 '4 MeV @. 

The analysis of these new data in kinematic conditions unaccessible in previous experiments, where Q 2 was less 
than 0.4 (GeV/c) 2 and T p generally around 100 MeV, requires a theoretical treatment where all relativistic effects are 
carefully and consistently included. Indeed these recent data are well described by RDWIA calculations JlJJll] |l8[ . 

Fully relativistic models based on the RDWIA have been developed by different groups [§-{L4j. It was shown in 
Ref. |12(| that the hadronic part of the relativistic transition amplitude can be written in terms of Schrodinger-like 
wave functions for bound and scattering states and of an effective nuclear current operator which contains the Dirac 
potentials. In this way, complications due to the solution of the Dirac equation can be avoided and spinor distortion can 
be accounted for solving a relativized Schrodinger equation, of the same type as that solved in ordinary nonrelativistic 
DWIA calculations with nonrelativistic or relativistic equivalent potentials, but Dirac scalar and vector potentials 
appear in the nuclear current operator. This so-called effective Pauli reduction does not represent an approximation. 
It is in principle exact and can be considered an alternative fully relativistic approach. 

This approach was adopted in Refs. [ |l9|]2"ofl with an effective momentum approximation (EM A) to incorporate 
spinor distortion into the effective current operator. In this approximation the momentum operators in the term a p, 
appearing in the lower components of the Dirac spinor and in the effective current operator, are replaced by their 
asymptotic values. This model contains approximations, but accounts for all the main relativistic effects and is able 
to give a good description of the most recent experimental results. The spectroscopic factor extracted in comparison 
with 16 0(e, e'p) data for lp-shell proton knockout is the same as in other RDWIA analyses, i.e. 0.7 Jl6[ . 

Both nonrelativistic DWIA and RDWIA models are able to describe, with a good degree of accuracy, the (e, e'p) data 
at low energies, but the nonrelativistic DWIA approach is more flexible. Some relativistic corrections can be included 
in a nonrelativistic treatment, but a fully relativistic model is needed for the analysis of the data at higher energies 
that are now becoming available. It is thus important to establish a clear relationship between the nonrelativistic 
DWIA treatment that was extensively used in the analysis of low-energy data and RDWIA treatments, in order to 
understand the relevance of genuine relativistic effects and the limit of validity of the nonrelativstic model. 

Relativistic effects as well as the differences between relativistic and norclativistic calculations have been widely and 
carefully investigated in different papers where RDWIA treatments have been developed. The differences, however, 
were usually investigated starting from the basis of a relativistic model where terms corresponding to relativistic 
effects are cancelled, such as, for instance, the lower components in the Dirac spinor and the Darwin factor, or where 
nonrelativistic approximations are included. Although very interesting, these investigations do not correspond to the 
result of a direct and consistent comparison between RDWIA and the DWIA calculations carried out with the program 
DWEEPY, that was used in the analysis of low-energy data. In fact, DWEEPY is based on a nonrelativistic treatment 
where some relativistic corrections are introduced: a relativistic kinematics is adopted and relativistic corrections at 
the lowest order in the inverse nucleon mass are included in the nuclear current operator ]2l| , p2| ], which is derived 
from the Foldy-Wouthuysen transformation [p3[ . The Darwin factor can be simply included in a nonrelativistic 
treatment |l5| , but it was not explicitly considered in the data analyses carried out with DWEEPY. 

Only indirect comparison between relativistic and norelativistic calculations can be obtained from the available data 
analyses carried out with DWEEPY and in RDWIA. In fact, the two types of calculations make generally use of different 
optical potentials and bound state wave functions, and the difference due to the different theoretical ingredients cannot 
be attributed to relativity. This problem was already clear in previous investigations of relativistic effects. Various 
attempts were made to reproduce the conditions of a nonrelativistic calculation from a relativistic approach, but a 
consistent and direct comparison has not yet been achieved. 

The main aim of this paper is to make clear the relationship between the DWIA approach in DWEEPY and RDWIA 
treatments, and to investigate the relevance of genuine relativistic effects through a direct comparison between the 
results of the two calculations. A fully relativistic RDWIA treatment of the (e, e'p) knockout reaction has thus been 
developed. The effective Pauli reduction [jl2| has been adopted for the scattering state. For the bound state the 
numerical solution of the Dirac equation has been used, as in this case it does not represent a too difficult problem. 
Various computer programs are in fact available, able to explain the global and s.p. properties of a nucleus within the 
framework of a relativistic mean-field theory. For the scattering state, where a partial-wave expansion is performed, 
the effective Pauli reduction appears simpler and more flexible, and it is equivalent to the solution of the Dirac 
equation. From this point of view, our relativistic approach does not contain approximations, since the Schrodinger- 
like equation is solved for each partial wave starting from a relativistic optical potential and without assuming the 
EMA prescription in the effective current operator. 

The numerical results of this RDWIA model have been compared with the results of DWEEPY. A consistent comparison 
requires the same kinematical conditions and the use of consistent theoretical ingredients in the two calculations. 
Thus, in DWEEPY we have adopted for the bound state the normalized upper component of the Dirac spinor and for 
the scattering state the solution of the same Schrodinger-equivalent optical potential of the relativistic calculation. 
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Different kinematics have been considered for the comparison, with the aim to investigate the relevance of relativistic 
effects not included in DWEEPY in different situations, at low energy, in the kinematical region where DWEEPY was 
extensively and successfully applied, and at higher energy, where relativistic effects are more evident, in order to 
establish the limit of validity of the nonrelativistic approach. 

The relevant formalism is outlined in Sec. II. Relativistic and nonrelativistic calculations of the cross section and 
response functions for the 16 0(e,e'p) reaction are compared in Sec. Ill, where various relativistic effects are investi- 
gated. Even though the comparison with experimental data is not the main aim of this work, in Sec. IV we check the 
reliability of our approach in comparison with data. Some conclusions are drawn in Sec. V. 



II. FORMALISM 



A. Relativistic current 



In the one-photon exchange approximation, where a photon is exchanged between the incident electron and the 
target nucleus, the coincidence cross section of the (e, e'p) reaction is given by the contraction between the lepton 
tensor, dependent only on the electron variables and completely determined by QED, and the hadron tensor, whose 
components are given by bilinear products of the transition matrix elements of the nuclear current operator. According 
to the impulse approximation, in which only one nucleon of the target is involved in the reaction, the nuclear current 
is assumed to be a one-body operator. 

In RDWIA the matrix elements of the nuclear current operator, i.e. 

J' 1 = J dr W/(r) ~f exp{iq • r} *i(r), (1) 

are calculated with relativistic wave functions for initial bound and final scattering states. 

We choose the electromagnetic current operator corresponding to the cc2 definition of Ref. [Q, i.e. 

T = fi(QV + i ^MQ 2 ) <r» v qu, (2) 

where q" — (u, q) is the four-momentum transfer, Q 2 = q 2 — uj 2 , F\ and F 2 are Dirac and Pauli nucleon form factors, k 
is the anomalous part of the magnetic moment, and <r Miy = i/2 [7^, 7"] . Current conservation is restored by replacing 
the longitudinal current by 

J L = J z = ~ J°. (3) 

In our reference frame the z-axis is along q and the y-axis parallel to q x p'. 
The bound state wave function 

•-(S) (4) 

is obtained as the Dirac-Hartree solution from a relativistic Lagrangian with scalar and vector potentials. Several 
computer codes calculating the ground and excited state properties of nuclei are easily available in literature |^5|,^6| . 

The ejectile wave function is obtained following the direct Pauli reduction method. It is well known that a Dirac 
spinor 

*=(*-) < 5 > 

can be written in terms of its positive energy component \&+ as 

V E+M+S-V + / 

where S = S(r) and V = V(r) are the scalar and vector potentials for the nucleon with energy E. The upper 
component can be related to a Schrodinger-like wave function $ by the Darwin factor D(r), i.e. 

vf + = V^w $, (7) 
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D(r) = 



E + M + S(r) -V(r) 
E + M ' 



(8) 



The two-component wave function $ is solution of a Schrodinger equation containing equivalent central and spin- 
orbit potentials, which are functions of the scalar and vector potentials S and V and are energy dependent. 

Inserting Eq. (|7]) into Eq. (j^) and using the relativistic normalization, we obtain the wave function for the knocked 
out nucleoli 



T xt o E' + M 
t f ' V 2E' 
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E' + M 



2E' 
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(9) 



where 



C(r) =E' + M + S(r)-V(r). 
If we substitute Eqs. (§), (§), and (|) into Eq. (0), we obtain the relativistic nuclear current 



(10) 
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exp{iq ■ r}, 



(11) 



where the operator p has been replaced by the gradient — iV, which operates not only on the components of the Dirac 
spinor but also on exp{iq • r}. 



B. Nonrelativistic current 



In nonrelativistic DWIA the matrix elements of the nuclear current are 

J£ r = J dr (r) f nr exp{iq • r} ^(r), (12) 

where the bound and scattering states are eigenfunctions of a Schrodinger equation. 

In standard DWIA analyses phenomenological ingredients are usually adopted for <5>i(r) and $/(r). In this work 
and in order to perform a consistent comparison with RDWIA calculations, we employ for $j(r) the upper component 
of the Dirac wave function Ui and for $/(r) the Schrodinger-like wave function that appear in the relativistic current 
ofEq. (U). 

The nuclear current operator is obtained from the Foldy-Wouthuysen reduction of the free-nucleon Dirac current 
through an expansion in a power series of l/M, i.e. 

&- = E3w- (13) 

In the program DWEEPY the expansion is truncated at second order (N = 2). 
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C. Cross section and response functions 



The coincidence cross section of the unpolarized (e, e'p) reaction is written in terms of four nuclear structure 
functions f\\i [|| as 



cr = cm E'\p'\ \ poo/oo + P11/11 + Poi/oi cosa + cos 2a }, 



(14) 



where ctm is the Mott cross section and a the out of plane angle between the electron-scattering plane and the (q,p')- 
plane. The coefficients p\\> are obtained from the components of the lepton tensor and depend only on the electron 
kinematics. The structure functions fw are given by suitable combinations of the components of the nuclear current 

as 



/ 00 = < J°(J ) f >, 

f n = < j x (ry > + < jy(jy) ] >, 

/oi = -2V2 ate r< j x (j°y >i, 
h-x = <j y {j y y >-<j x (j x Y >, 



(15) 



where the brackets mean that the matrix elements are averaged over the initial and summed over the final states 
fulfilling energy conservation. 

In the following we use a different definition of the structure functions, that is 



Rl = (2^) 3 /oo, 

RhT ~ ~7T /o1 ' 



R T = (2^) 3 /n, 
Rtt = -(27T) 3 /l-l- 



(16) 



If the electron beam is longitudinally polarized with helicity h, the coincidence cross section for a knocked out 
nucleon with spin directed along s can be written as l| 



1 



1 + P-s + h 



[A + P' ■ s) 



(17) 



where Co is the unpolarized cross section of Eq. ( |l4| ) , P the induced polarization, A the electron analyzing power and 
P' the polarization transfer coefficient. We choose for the polarimeter the three perpendicular directions: L parallel 
to p', N along qxp', and T = N x L. The corresponding polarization observables can be written in terms of new 
structure functions, which contain explicitly the polarization direction of the emitted nucleon. In coplanar kinematics 
(a = 0,7r), only P N , P' L , and P' T survive §. 



III. RELATIVISTIC EFFECTS IN THE 16 0(e,e'p) REACTION 

The reaction 16 0(e, e'p) 15 N has been chosen as a well suited process for testing the relativistic program and 
investigating the differences with respect to the nonrelativistic program DWEEPY H. A large experience concerning 
theoretical calculations is available on this reaction and a considerable amount of experimental data at different 
energies and kinematics has been published, including polarization measurements. 

The relativistic bound-state wave functions used in the calculations have been obtained from the program ADFX 
of Ref. p6j, where relativistic Hartree-Bogoliubov equations are solved. The model is applied in the mean-field 
approximation to the description of ground-state properties of spherical nuclei f27j] . Sigma-meson, omega- meson, 
rho-meson and photon field contribute to the interaction and their potentials are obtained by solving self-consistently 
Klein-Gordon equations. Moreover, finite range interactions are included to describe pairing correlations and the 
coupling to particle continuum states. The corresponding wave function for the nonrelativistic calculation has been 
taken as the upper component of the relativistic Dirac four-component spinor with the proper normalization, i.e. 
normalized to one in the coordinate and spin space. This is not the best choice for DWEEPY, but a consistent use of 
the theoretical ingredients is necessary to allow a clear comparison between the two approaches. 

The relativistic nuclear current was taken as in Eq. (^) p4| . This expression is not only more fundamental than the 
other forms recovered from the Gordon decomposition, but it is also consistent with the nonrelativistic current used 
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in DWEEPY, where the relativistic corrections up to order 1/M 2 are obtained from a Foldy-Wouthuysen transformation 
applied to the interaction Hamiltonian where the nuclear current has the same form as in Eq. (^|). The Dirac and 
Pauli form factors F\ and F2 are taken in both calculations from Ref. p8| ]. 

The outgoing-proton wave function is calculated by means of the relativistic energy-dependent optical potential of 
Rcf. |2jJ, which fits proton elastic scattering data on several nuclei in an energy range up to 1040 MeV. The program 
GLOBAL of Ref. > which generates the scalar and vector components of the Dirac optical potential, has been used. 
Different fits, available from the code, were explored. The Schrodinger equivalent potentials calculated in the same 
program were used for the nonrelativistic calculation. 

The comparison between the results of the two approaches is not restricted to the cross section, but involves also 
the structure functions, which can be experimentally separated and show a different sensitivity to the treatment of 
the theoretical ingredients and to relativistic effects. 

A kinematics with constant (q,w), or perpendicular kinematics, was chosen as convenient for the comparison, but 
some calculations were performed also in parallel kinematics. In both kinematics the incident electron energy and the 
outgoing proton energy are fixed. In the kinematics with constant (q,u)) the electron scattering angle is calculated by 
imposing the condition |q| = \p'\. Changing the angle of the outgoing proton, different values of the recoil or missing 
momentum p m , with p nl = q — p', can be explored. In parallel kinematics q is parallel or antiparallel to p' , and 
different values of p m are obtained changing the electron scattering angle, and thus q. 

The beam energy Eq was fixed in the present work at 2 GeV, in order to minimize the effect of the Coulomb 
distortion, which is included in the relativistic program only through the effective momentum approximation [^,0. 
Calculations were performed for an outgoing proton energy up to 400 Mev. The missing momentum was explored in 
a range up to ^400 MeV/c. 



A. Relativistic vs nonrelativistic results 



In this Section the results of the comparison between the DWIA calculations performed with DWEEPY and RDWIA 
calculations are discussed. DWEEPY is based on a nonrelativistic treatment, but does already contain some relativistic 
corrections in the kinematics and in the nuclear current through the expansion in 1/M. Therefore, the results of 
DWEEPY cannot be obtained from the relativistic program simply by dropping relativistic effects, such as the lower 
components of the Dirac spinor and applying the proper normalizations. Here the comparison is done between the 
results of the two independent programs. In the first place we checked the numerical consistency of the two programs 
and verified that they give the same result in the same situation, i.e. when all the differences are eliminated. This 
numerical check gives us confidence that we are not interpreting the contribution of different ingredients as a relativistic 
effect. 

The comparison between the results of the new RDWIA program and DWEEPY is shown in Figs. 1-6, for the structure 
functions and the cross section of the 16 0(e, e'p) 15 N g s . reaction in a kinematics with constant (q,u>), at three values 
of the proton energy, i.e. T p = 100 MeV, 200 MeV and 300 MeV. 

It is clear from these figures that the differences rapidly increase with the energy. Moreover, the relativistic result 
is smaller than the nonrelativistic one. This is a well known effect, which was found in all the relativistic calculations 
and which is essentially due to the Darwin factor jljj of Eq. (^). In addition, the relativistic calculations include the 
typical normalization factor (E + M)/2E (see Eq. which has the value 0.95 at 100 MeV and decreases to 0.87 
at 300 MeV. 

A significant difference is found for the transverse structure function i?x even at T p = 100 MeV, where a reduction 
of about 15% is obtained with respect to the nonrelativistic calculation. The reduction grows up to about 25% at 200 
MeV, and 40% at 300 MeV. The difference is sensibly reduced, mainly at lower energies, by including in the nuclear 
current the terms to the order 1/M 3 . 

Only small differences arc found for the longitudinal structure function i?L at all the considered proton energies. 
Its size, however, decreases when the energy increases and therefore its contribution to the cross section becomes less 
important. 

Large differences are generally found for the interference structure function Rut- The combined relativistic effects 
on i?x and Rjjt are responsible for the different asymmetry in the cross section at positive and negative values of p m , 
where the mismatch between the two structure functions adds up on the one side and compensates on the other side. 

Large relativistic effects are also found on Rtt, which is anyhow much smaller than the other structure functions 
and gives only a negligible contribution to the cross section. 

As relativistic effects increase with the energy, the conclusion of this comparison is that DWEEPY can be used with 
enough confidence at energies around 100 MeV, and, with some caution, up to about 200 MeV. Higher-order terms in 
the nuclear current can account for a part of the difference in this kincmatical region. A fully relativistic calculation is 
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anyhow convenient at 200 MeV and necessary above 300 MeV. This result is consistent with the old one of Ref. p2] , 
where the validity of the relativistic corrections to the nuclear current, calculated as an expansion on 1/M, were 
discussed and an upper limit of |g| ~ 600 MeV/c was stated for the nonrelativistic calculations. 

A calculation was performed also in parallel kinematics at T p — 100 and 200 MeV. In this kinematics only the 
longitudinal i?L and the transverse structure function i?x survive. The results are shown in Figs. 7 and 8. Also in 
this case small differences are found for i?L, while Rt is significantly reduced in RDWIA. The reduction increases 
with the outgoing proton energy, and, at a given values of T p , is larger at positive p m , when the momentum transfer 
decreases. This is mainly due to the fact that the leading term of the spin current, proportional to q, is the same in 
the relativistic and nonrelativistic expressions. 



B. Darwin factor and spinor distortion of the scattering state 

In this Section we shall discuss the effect of the optical potential in the Pauli reduction of the four-component Dirac 
spinor for the scattering state. We do not discuss the corresponding effect on the bound state, as in our calculations 
the bound-state wave function is taken directly from the solution of the relativistic Dirac equation. 

The effect of the optical potential on the Pauli reduction is twofold: the Darwin factor D(r) of Eq. (||), which 
directly multiplies the Schrodinger-equivalent eigenfunction, and the spinor distortion C (r) of Eq. (|l^) , which applies 
only to the lower component of the Dirac spinor. The distortion of the scattering wavefunction, which is calculated 
through a partial wave expansion, is always included in the calculation and affects in a similar way the relativistic 
and nonrelativistic result. 

The Darwin factor gives a reduction of about 5-10% at 100 MeV, which is consistent with the qualitative prediction 
of Ref. jyl . On the contrary, spinor distortion produces an enhancement of the cross section, so that the combined 
effect of the two corrections is in general small. Moreover, it decreases with the energy. The result is qualitatively in 
agreement with the ones of Ref. 0,0 ■ 

In Figs. 9 and 10 the comparison between the two results is shown in the kinematics with constant (q,u>) at 
T p = 100 MeV, for the structure functions and the cross section, respectively. The effect is always small. Spinor 
distortion enhances the cross section at high and negative recoil momenta, but this effect seems absent at p m > or 
it is pushed to very high momenta. 

The effect is larger on R^ than on Rt- Thus, at higher energies, where Rt becomes larger than i?L, the effect of 
spinor distortion on the cross section decreases. 



C. Effective momentum approximation 

In this Section the validity of EMA is discussed. This prescription, which consists in evaluating the momentum 
operator in the effective nuclear current using the asymptotic value of the outgoing proton momentum, simplifies 



considerably the numerical calculations, avoiding the evaluation of the gradient in Eq. (11). It is exact in plane wave 
impulse approximation (PWIA), where the scattering wave functions are eigenfunctions of the momentum, but in 
DWIA it disregards the spreading of the distorted proton wave function in momentum space due to FSI. 

This approximation was used in some relativistic calculations, and in particular in the model of Refs. []l9| , po[ for 
bound and scattering states. Since in our approach the bound-state wave function is taken directly from the solution 
of the Dirac equation, we investigate the validity of EMA only for the scattering state. 

We have to notice that in our calculations EMA does not change the nuclear current operator, which is calculated 
with the cc2 formula and therefore does not depend on the momentum. The only dependence is contained in the 
Pauli reduction of the scattering wave function. 

The effect of EMA in our RDWIA approach is shown in Figs. 9 and 10, where the the structure functions and 
the cross section calculated with EMA in the kinematics with constant (q, w) at T p = 100 MeV are displayed and 
compared with the exact result. At 100 MeV the difference is indeed sensible, but it rapidly decreases with the 
energy. At 200 MeV it is much smaller, and becomes really negligible at ~ 400 MeV, and therefore at the energy of 
the TJNAF experiment. This behaviour can be understood if one considers that distortion effects decrease with the 
energy, so that at high energy DWIA approaches the PWIA result, where EMA is exact. 
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IV. COMPARISON WITH EXPERIMENTAL DATA 



In this Section we shall discuss the comparison of our RWDIA results with experimental data. Data are available 
at low energies, where DWEEPY was extensively and successfully applied, and, more recently, at higher energies, where 
other RDWIA calculations have given an excellent agreement. Even though a precise description of experimental data 
is not the main aim of this work, it is interesting to test the predictions of our relativistic approach in comparison 
with data that have been already and successfully described by other theoretical models. 



A. The 1(i O(e,e'p) reaction 

Low-energy data are presented in terms of the reduced cross section || , which is defined as the cross section divided 
by a kinematical factor and the elementary off-shell electron-proton scattering cross section, usually a cc \ of Ref . [^4) . 

In Fig. [ll] the reduced cross sections measured at NIKHEF |3(J for the le O(e,e'p) knockout reaction and for the 
transitions to the l/2~ ground state and the 3/2~ excited state of 15 N are displayed and compared with the results 
given by our RDWIA program and by DWEEPY. The experiment was carried out in parallel kinematics at T p = 90 
MeV. 

The relativistic results are lower than the nonrelativistic ones and the corresponding spectroscopic factors are 
approximately 10% larger than those deduced from nonrelativistic analyses. Thus, the normalization (spectroscopic) 
factor, applied in Fig. [n] to the calculated results in order to reproduce the size of the experimental data, is 0.70 for 
RDWIA and 0.65 for DWEEPY. The same value is adopted for the two final states. 



As we already stated in Sec. HI, only small differences are found at this proton energy between the two models. 
Thus, the results of the two calculations are almost equivalent in comparison with data, which are reasonably described 
by both calculations. A better agreement is found for the l/2~ than for the 3/2~ state. In any case, it is not as 
good as in the DWIA analysis of Ref. J3(j performed with DWEEPY. This result is expected, as we already claimed that 
the theoretical ingredients, bound state and optical potential, used in the present calculation do not represent the 
best choice for DWEEPY, but are here adopted in order to allow a consistent comparison between the relativistic and 
nonrclativtsic approaches. 

In Figs. |||-|l5|the cross sections and the structure functions measured at TJNAF [|l6) for the 16 0(e,e'p) knockout 
reaction and for the transitions to the l/2~ ground state and the 3/2~ excited state of 15 N are displayed and compared 
with the results of our RDWIA model. The experiment was carried out in a kinematics with constant (q, lu), with Eq = 
2.4 GeV and cu ~ 439 MeV. Only RDWIA calculations are shown in the figure, since we know from the investigation of 



Sec. Ill that at the proton energy of this experiment relativistic effects are large and a relativistic analysis is necessary. 
In order to study the sensitivity to different optical potentials, we compare in the figures results obtained with the 
EDAD1 and EDAI-0 fits of Ref. |2S|j . Only small differences are given by the two optical potentials. The agreement 
with data is satisfactory and of about the same quality as in other RDWIA analyses, but for the interference structure 
function i?LT for the l/2~ state at intermediate missing momenta. A better description of data might be obtained 
with a more careful determination of the theoretical ingredients. 

It is interesting to notice that the spectroscopic factor applied to all the calculatations in Figs. |l2|-|l5| is the same, 
i.e. 0.7, as that found in the comparison with NIKHEF data shown in Fi g. O . A spectroscopic factor of about 0.7 
was also obtained in previous RDWIA analyses of the same TJNAF data |16|. 



B. The 12 C(e,e'p) and le O(e,e'p) reactions 

The induced polarization of the outgoing proton P N was measured at Bates for the 12 C(e, e'p) reaction [Q. Data 
were taken in a kinematics with constant (q, lo) at Eq = 579 MeV and to ~ 290 MeV. In Fig. these data are displayed 
and compared with our RDWIA calculations. Results obtained with the EDAD1 and EDAI-C optical potentials are 
compared in the figure. The EDAD1 curve gives a better description of the experimental data at high values of the 
missing momentum, but both calculations are in fair agreement with data. With EDAD1 potential, we also plot 
results after eliminating the negative energy components in the bound state. 

As already shown in Ref. [^lj the polarization P N is enhanced by the presence of the negative energy components 
of the relativistic bound state wave function. This result is confirmed in Fig. [l6| by the dashed line which gives a 
smaller P N . 

A slightly higher polarization is obtained in Ref. with the nuclear current written in the ecl form according to 
Ref. |§. 
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The components of the polarization coefficient P' L and P' s were measured at TJNAF [jl7j for the 16 0(e, e'p) reaction 
and for the transitions to the ground state and the 3/2~ excited state of 15 N. The experiment was performed 

], that is the one of Figs. [l2|-p^[ The experimental data are compared with our 
RDWIA calculations in Fig. |T7j for l/2~and in Fig. |l| for the 3/2" state. The two curves in the figures show the 
results obtained with the EDAD1 and EDAFO fits. Both results are in satisfactory agreement with data. 



in the same kinematics as in Rcf. 16 



V. SUMMARY AND CONCLUSIONS 



In order to clarify the differences between the usual nonrelativistic approach and a relativistic description of exclusive 
(e, e'p) knockout reactions, we have performed a fully relativistic calculation and compared it with the nonrelativistic 
results of the DWEEPY code, that was successfully used to analyse a large number of experimental data. The transition 
matrix element of the nuclear current operator is written in RDWIA using the relativistic bound state wave functions 
obtained in the framework of the mean field theory, and the direct Pauli reduction method with scalar and vector 
potentials for the ejectile wave functions. Correspondingly, the nonrelativistic DWIA matrix elements are computed 
in a consistent way to allow a direct comparison with the relativistic results. 

The main aim of this paper was not to make a precise analysis of the existing experimental data, but to discuss 
the use of a nonrelativistic approach at energies higher than those generally considered up to few years ago, and to 
clarify the possible relativistic effects arising also at lower energies. 

We have used the new RDWIA and DWEEPY codes to perform calculations for several kinematical conditions. The 
relativistic results are always smaller than the nonrelativistic ones and the difference increases with energy. The 
transverse and interference structure functions are particularly sensitive to relativistic effects, much more than the 
longitudinal structure function. Rt is sensibly reduced even at low energy: inclusion of higher order terms in the 
nonrelativistic nuclear current can reduce the difference, but a fully relativistic calculation is necessary above T p ~ 
200 MeV. 

The effect of the scalar and vector potentials in the Pauli reduction for the scattering state has been discussed. These 
potentials appear in the relativistic treatment and are absent in the nonrelativistic one. The combined contribution 
of the reduction due to the Darwin factor and of spinor distortion, which enhances the effects of the lower components 
of the Dirac spinor, is always small. 

The validity of EMA in the scattering state of relativistic calculations has been studied. The differences with respect 
to the exact results are sensible at T p = 100 Mev, but rapidly decrease with the energy and become negligible at T p ~ 
400 MeV. 

We have tested our new RDWIA calculations in comparison with experimental data that have already been described 
by other models. The agreement is satisfactory and comparable with other relativistic analyses. 



[1] S. Frullani and J. Mougey, Adv. Nucl. Phys. 14, 1 (1984). 

[2] S. Boffi, C. Giusti, and F. D. Pacati, Phys. Rep. 226, 1 (1993) 

[3] S. Boffi, C. Giusti, F. D. Pacati, and M. Radici, Electromagnetic Response of Atomic Nuclei, Oxford Studies in Nuclear 

Physics (Clarendon Press, Oxford, 1996). 
[4] J. J. Kelly, Adv. Nucl. Phys. 23, 75 (1996). 
[5] J. Mougey et al. Nucl. Phys. A262, 461 (1976); 

M. Bernheim et al. Nucl. Phys. A375, 381 (1982). 
[6] P. K. A. de Witt Huberts, J. Phys. G 16, 507 (1990); 

L. Lapikas, Nucl. Phys. A553, 297c (1993). 
[7] C. Giusti and F. D. Pacati, Nucl. Phys. A473, 717 (1987); Nucl. Phys. A485, 461 (1988). 
[8] A. Picklesimer, J. .W .Van Orden, and S. J. Wallace, Phys. Rev. C 32, 1312 (1985); 

A. Picklesimer and J. .W Van Orden, Phys. Rev. C 40, 290 (1989). 
[9] Y. Jin, D. S. Onley, and L. E. Wright, Phys. Rev. C 45, 1311 (1992); 

Y. Jin, J. K. Zhang, D. S. Onley, and L. E. Wright, Phys. Rev. C 47, 2024 (1993); 

Y. Jin and D. S. Onley, Phys. Rev. C 50, 377 (1994). 
[10] J. M. Udfas, P. Sarriguren, E. Moya de Guerra, E. Garrido, and J. A. Caballero, Phys. Rev. C 48, 2731 (1993); 

J. M. Udfas, P. Sarriguren, E. Moya de Guerra, and J. A. Caballero, Phys. Rev. C 53, R1488 (1996); 

J. M. Udfas, J. A. Caballero, E. Moya de Guerra, J. E. Amaro, and T. W. Donnelly, Phys. Rev. Lett. 83, 5451 (1999). 







J. M. Udfas and R. Vignote Phys. Rev. C 62, 034302 (2000). 

M. Hedayati-Poor, J. I. Johansson, and H. S. Sherif, Phys. Rev. C 51, 2044 (1995). 
J. I. Johansson, H. S. Sherif, and G. M. Lotz, Nucl. Phys. A605, 517 (1996). 
J. I. Johansson and H. S. Sherif, Phys. Rev. C 59, 3481 (1999). 

S. Boffi, C. Giusti, F. D. Pacati, and F. Cannata, Nuovo Cimento A 98, 291 (1987). 

J. Gao et al. Phys. Rev. Lett. 84, 3265 (2000). 

S. Malov et al. Phys. Rev. C 62, 057302 (2000). 

R. J. Woo et al. Phys. Rev. Lett. 80, 456 (1998). 

J. J. Kelly, Phys. Rev. C 56, 2672 (1997); 59, 3256 (1999). 

J. J. Kelly, Phys. Rev. C 60, 044609 (1999). 

K. W. McVoy and L. Van Hove, Phys. Rev. 125, 1034 (1962). 

C. Giusti and F. D. Pacati, Nucl. Phys. A336, 427 (1980). 

L. L. Foldy and S. A. Wouthuysen, Phys. Rev. 78, 29 (1950). 

T. de Forest, Jr., Nucl. Phys. A392, 232 (1983). 

C. J. Horowitz, D. P. Murdoch, and B. D. Serot, in Computational Nuclear Physics I: Nuclear Structure, edited by 

K. Langanke, J. A. Maruhn, and S. E. Koonin (Springer- Verlag, Berlin, 1991). 

W. Poschl, D. Vretenar, and P. Ring, Comput. Phys. Commun. 103, 217 (1997). 

G. A. Lalazissis, J. K6nig, and P. Ring, Phys. Rev. C 55, 540 (1997). 

P. Mergell, Ulf-G. Meissner, and D. Drechsel, Nucl. Phys. A596, 367 (1996). 

E. D. Cooper, S. Hama, B. C. Clark, and R. L. Mercer, Phys. Rev. C 47, 297 (1993). 

M. Leuschner et al, Phys. Rev. C 49, 955 (1994). 



10 




100 200 300 400 100 200 300 400 

p m [MeV/cl p m [MeV/c] 



FIG. 1. The structure functions of the ie O(e,e'p) reaction as a function of the recoil momentum p m for the transition to 
the 1/2" ground state of 15 N, in a kinematics with constant (q,uj), with Eq — 2000 MeV and T p = 100 MeV. The solid lines 
give the RDWIA result, the dotted lines the nonrelativistic result of DWEEPY. Positive (negative) values of p m refer to situations 
where the angle between p' and the incident electron k is larger (smaller) than the angle between q and k. 
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FIG. 2. The same as in Fig. 1, but for the cross section of the le O(e,e'p) reaction. 
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FIG. 3. The same as in Fig. 1 at T p = 200 MeV. 
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FIG. 4. The same as in Fig. 2 at T p = 200 MeV. 
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FIG. 5. The same as in Fig. 1 at T p = 300 MeV. 
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FIG. 6. The same as in Fig. 2 at T p = 300 MeV. 
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FIG. 7. The longitudinal and transverse structure functions of the ie O(e, e'p) reaction as a function of the recoil momentum 
p m for the transition to the l/2~ ground state of 15 N, in parallel kinematics with Eo = 2000 MeV and T p = 100 MeV. Positive 
(negative) values of jj m refer to situations where \p'\ > \q\ (\p'\ < \q\). Line convention as in Fig. 1. 
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FIG. 8. The same as in Fig. 7 for T p = 200 MeV. 
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FIG. 9. The structure functions of the ie O(e,e'p) reaction as a function of the recoil momentum p m for the transition to 
the l/2~ ground state of 15 N, in the same situation as in Fig. 1. The solid lines give the RDWIA result, the dotted lines the 
calculation without the Darwin factor and spinor distortion and the dashed line the EMA. 





FIG. 11. The reduced cross section of the ie O(e,e'p) reaction as a function of the recoil momentum p m for the transitions 
to the l/2~ ground state and to the 3/2~ excited state of lj N, in parallel kinematics with Eo = 520 MeV and T p = 90 MeV. 
The data are from Ref. po| . The solid lines give the RDWIA result, the dotted lines the nonrelativistic result of DWEEPY. 
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FIG. 12. The cross section of the 16 0(e, e'p) reaction as a function of the recoil momentum p m for the transition to the l/2~ 
ground state of 15 N in a kinematics with constant (q, uS), with _Bo = 2445 MeV and T p — 433 MeV. The data are from Ref. 
Jlgl . The solid line gives the RDWIA result with the EDAD1 optical potential, the dotted line the RDWIA result with the 
EDAI-0 optical potential. 
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FIG. 13. The same as in Fig. [l^ but for the response functions of the 16 0(e,e'p) reaction. 
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FIG. 16. The induced polarization of the emitted proton for the 12 C(e,e'p) reaction as a function of the recoil momentum 
p m for the transition to the 3/2~ ground state of n B in a kinematics with constant (q,uj), with Eo = 579 MeV and T p = 274 
MeV. The data are from Ref. jlg|| . The solid line gives the RDWIA result with the EDAD1 optical potential, the dotted line 
the RDWIA result with the EDAI-C optical potential, and the dashed line the result with the EDAD1 optical potential and 
after removing the negative energy components in the bound state. 
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FIG. 17. The components of the polarization transfer coeffcient P' L and P' T for the ia O(e, e'p) reaction as a function of the 
recoil momentum p m for the transition to the l/2~ ground state of 15 N in the same kinematics as in Fig. [jjj The data are 
from Ref. The solid lines give the RDWIA result with the EDAD1 optical potential, the dotted lines the RDWIA result 
with the EDAI-O optical potential. 
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FIG. 18. The same as in Fig. [Tt] but for the transition to the 3/2 excited state of 15 N. 
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